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The historical interplay between societies are governed by many factors, including in particular 
spreading of languages, religion and other symbolic traits. Cultural development, in turn, is coupled 
to emergence and maintenance of information spreading. Strong centralized cultures exist thanks 
to attention from their members, which faithfulness in turn relies on supply of information. Here, 
we discuss a culture evolution model on a planar geometry that takes into account aspects of the 
feedback between information spreading and its maintenance. Features of model are highlighted by 
comparing it to cultural spreading in ancient and medieval Europe, where it in particular suggests 
that long lived centers should be located in geographically remote regions. 
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I. INTRODUCTION 

The expansion and decline of social structures depend 
on information spreading in form of languages, religion, 
or other cultural inventions [1] . In the recent years many 
mathematical models have been proposed for social inter- 
actions and dynamics, trying to understand social struc- 
tures. However, a main feature of most models of culture 
dissemination is an adaptation toward local or global con- 
sensus [U |2) , an equilibration which is also found in voter 
models 4J, social impact theory [5], majority rules [B], 
the Sznajds model [7], the Deffuant model [8] and the 
bounded confidence models [?J]. 

The multitude of models that emphasize consensus- 
dynamics contrast a reality where consensus is often bro- 
ken by emergence of new cultures, languages or opinions. 
One driving force for heterogeneity is the need for atten- 
tion, where individuals not only aim at mutual under- 
standing, but at the same time also fight for individual 
attention. This attention battle is more than random 
fluctuations of agents [TU] or rejection of other opinions 
[llj : The battle may, for example, involve positive feed- 
back mechanisms as suggested by [Hj. In this paper, we 
take the possibility for a new culture to emerge into ac- 
count, in addition to the local alignment rules. Lacking a 
simple realistic mechanism for creation of new cultures, 
we here simply parametrize this "emergence" in terms of 
a rate p n ew for initiation of new cultures. 

Another common theme of models dealing with the 
spread of information is that two different pieces of in- 
formation are treated on an equal footing. This is in 
general an incorrect assumption, as the importance of 
two bits of information in general is asymmetric. One 
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FIG. 1: (Color on-line) System snapshots at various times 
t = {47,48,49} x 10 3 (from left to right), i.e. consecutive 
columns present subsequent states of the system separated by 
1 000 Monte Carlo steps. Different panels demonstrates: spa- 
tial structure of cultures (top panel), distance to the cultural 
center (middle panel) and information pathways (bottom 
panel). Cultural centers are located on the "white- blue stars" 
in the middle panels. Parameters: p = 0.01, p ncw = 0.004 and 

Prcpcat = 0.015. 



sorting principle is to use the information age as a sort- 
ing criteria, reflecting the fact that the value of infor- 
mation typically decays with time |13j . Previous studies 
P31 HI] demonstrated that such a sorting principle has 
major consequences for the spatio-temporal dynamics of 
information. Importantly, newer information overriding 
older information has been observed in spreading of lin- 
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guistic features [151 HZ]- Furthermore, a simple model of 
diffusion of information from a cultural stronghold with 
age sorting is shown to be compatible with the observed 
pattern of word distribution in Japan |18| . 

When individuals sort information based on its age, an 
existing cultural center will continuously need to gener- 
ate new information to maintain their sphere of influence. 
Accordingly, we here characterize the strength of a cul- 
tural center by the rate with which it is able to generate 
fashions, Repeat- We will take this rate as a characteristic 
of a cultural center, and it keeps generating new fashions 
until the center is eliminated by information generated 
from competing cultural centers. 

A main feature of the famous Axelrods model [21 |3] 
for social alignment is conservativeness in communi- 
cation, implemented by having individuals with many 
types of opinions and a preference for communication 
between individuals that share many traits. This pref- 
erence makes people more open for communication to- 
wards sources where they earlier obtained information. 
We here parametrize such conservativeness into a single 
parameter, p, that counts the chance that a given "site" 
or "agent" changes who he prefers to obtain information 
from. 

Overall our model aim to discuss the information flow 
association to emerging and collapsing cultural centers, 
each influencing their surroundings by an ongoing gen- 
eration of announcements that maintain their sphere of 
influence. The details of the model are presented in the 
next section. 



II. MODEL 

Our model considers many rumors / fashions / view- 
points / stories / ideas (denoted fashions in the follow- 
ing) competing on a two dimensional square lattice of 
sites which at any given moment can be occupied by one 
fashion only. We imagine each site as an agent, which 
in fact could be a whole group of people that by defini- 
tion share the same taste. Each site listens to their im- 
mediate four neighboring sites with a history dependent 
frequency. When communicating, they accept fashions 
only when they are newer than the current local fashion. 
History, or conservativeness, is quantified in terms of a 
preference in listening towards the direction where the 
last new idea came from. This is parametrized by the 
probability p (p < 0.25) to listen to one of the other 3 
directions. 

The fashions in the system come from cultural centers. 
Contrary to our previous model where only one cultural 
stronghold is placed in the system a priori |18j . we here 
assume that a new cultural center can emerge at any site 
with a small probability p nC w Such a site is recognized 
as a cultural center as long as it has its own fashion and 
not invaded by fashions from other sites. An existing 
cultural center in addition broadcasts itself repeatedly 
by initiating a new fashion with a rate p rep eat ■ 



We perform Monte-Carlo simulation of the model with 
parallel (synchronous) update of all agents. In the model, 
a fashion /(s, a) is characterized by its center s (the site 
at which the idea started) and its age a (how long time 
ago the fashion originated at the center). At any time t 
each site i have its current fashion Ii{t) — I(si,a,i) and 
its preferred direction from which this fashion was 

obtained. 

The time step from time t to t + 1 consists of the fol- 
lowing procedures: 

(i) Emergence of new cultural centers. With a proba- 
bility p n ew , a site i is randomly chosen out of all N 
sites in the system to become a new cultural cen- 
ter. It starts its own new fashion, i.e., Ii(t) is set 
to l(i,0). 

(ii) Repeated broadcast by existing cultural centers. 
Each cultural center i (i.e., = i) will start 
to spread a new fashion with probability p rep eat, 
namely Ii(t) is set to be I(i, 0). Putting it dif- 
ferently, every cultural center can re-broadcast the 
same fashion as a new one making it more appeal- 
ing. 

(iii) Spreading of ideas. For each site i in the system, 
the preferred site di{t) is chosen with probability 
1 — 3p (p < 0.25), or alternately one chooses one of 
the other neighbor sites with the probability p. The 
age cik of fashion Ik(t) = I(sk, afc) at the chosen site 
fc is compared with the age of the fashion at the 
site i. If a*, < a,, the site i accepts the fashion from 
the site k, namely set Ii(t + 1) = Ik(t) = I(sk,cik) 
and update its preferred direction to di(t + 1) = k. 
Otherwise the site i keep its fashion unchanged, i.e. 
Ii(t + 1) = Ii{t) = I(si,a,i) and keep its preferred 
direction di(t+ 1) = di(t). If a site i was a cultural 
center, the acceptance of competing idea destroys 
its ability as a cultural center, hence it stops re- 
peatedly broadcasting new fashions. 

(iv) Update of time. The ages of all fashions on all sites 
are increased by one. 

We simulate the model on an L x L square lattice under 
the periodic boundary conditions. We will also consider 
the model on a map of Europe, where the closed bound- 
ary conditions are imposed toward the sea regions. Initial 
condition is set so that no one has opinion nor preferred 
directions (all the sites weighted equally). We only in- 
vestigate properties of the system after the number of 
cultural centers have reached the steady state value. 

III. RESULTS 

1. Dynamics 

Dynamics of the model is depicted in Fig. [T] which 
presents time evolution on the European map. In the 
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top panel the respective cultures are presented with dif- 
ferent colors. The middle panel show current distance 
of each sites from its respective cultural center, defined 
as the place where its current fashion was introduced. 
The consecutive images illustrate the dynamics of the 
system, with meanderings of borders, as well as emer- 
gence of new cultural centers and disappearance of oth- 
ers. Bottom panel presents the information pathways 
(river landscape), based on the preferred direction for 
each site di(t). These arrows define the path to the cul- 
tural center to which every agent belongs. Intensity of 
points in the river landscape indicates how many times 
information was transmitted through every node, i.e. ev- 
ery time the idea is copied the number of transitions on 
all the nodes on the path (up to the origin) are increased 
by one. It shows clear river basin structure centered 
around respective cultural centers, much like what was 
obtained for the word spreading model of Ref. [18] . 
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FIG. 2: (Color on-line) Phase diagram showing average age, 
(age), (left column) and average size, (size), (right column) 
of cultural centers as the function of the probability of intro- 
duction of new cultures p ncw and the probability of resending 
new signals/rumor p rcp oat. The system size is 100 x 100. The 
top row corresponds to p = 0.25, while the bottom presents 
results for p = 0.01. 



2. Analysis of parameters 

The role of the three parameters p, p n0 w, and p r0 pcat is 
summarized in Fig. [2] representing behavior of 100 x 100 
system with periodic boundary conditions. In Fig. [2] the 
average age, (age), (left panel) and average size, (size), 
(right panel) of cultural centers are presented as a func- 
tion of Pncw and Prepeat, with p — 0.25 (top panel) and 
p = 0.01 (bottom panel) respectively. Note that p = 0.25 
examines the case where there is no conservativeness in 



the dynamics. As p becomes smaller, direction to exist- 
ing cultural centers are preferred and it becomes difficult 
for a new cultural center to emerge, as can be seen in the 
longer lived cultural centers for smaller p. 

For a given p, cultures live longer for smaller p ne w and 
larger p r0 pcat, because small p ne w decreases emergence of 
new cultural centers, while large p rep eat assures stability 
of the cultural center, i.e. constant broadcasting of new 
signals is vital for maintenance of cultural centers. In the 
opposite limit of large p ncw and small p repe at the ongoing 
strong competition between various ideas results in short 
living cultural centers. 

The biggest cultural centers are observed for small 
Pnew and small p rcp eat, see the right column of Fig. [2] 
Such a combination of parameters reduces the compe- 
tition among cultural centers and guarantees that each 
culture has enough time to spread over the whole sys- 
tem, resulting in one dominating culture of the system 
size. The opposite limit of large p n c W and large p rep eat, 
on the other hand, means frequent emergence of new cul- 
tural centers which survives relatively well, resulting in 
the coexistence of many small cultural centers. 

We checked the effect of the system size by comparing 
the results with 200 x 200 system, and confirmed that 
the data in Fig. [2] collapses onto the data from the big- 
ger system very well as long as the data with the same 
P, Prepeat > and p ne w / N are compared. The only exceptions 
are observed when the average size of fashion reaches the 
system size (data not shown). 



3. Competition between cultural centers 

The interesting aspect of our model is the ongoing re- 
placement of old cultural centers with new ones, a dy- 
namics primarily governed by the emergence of new cul- 
tural centers. To quantify this, we examine where new 
cultural centers tend to emerge, when there are already 
established cultural centers in the system. We re-run 
history multiple times with using a given snapshot of a 
100 x 100 system as an initial configuration. For this ini- 
tial configuration, we insert a new cultural center to a 
site at {x, y) at time zero, and run the simulation to see 
how long the new center survives, i.e. until it is overwrit- 
ten by fashions from neighbor centers. This procedure is 
repeated 25 times for each site to estimate the survival 
probability as a function of time t. 

The top panel of Fig. [3] presents the 2D map of the 
survival probability after 500 time steps (left map) and 
10 000 time steps (right panel) as a function of the po- 
sition (x,y) where new cultural center was inserted. In 
the plot, the positions of the cultural centers that exist 
in the initial configuration are marked by arrows. We 
see that new cultural centers are successful either when 
emerging very close to the existing cultural centers or 
when exploring remote regions. 

In order to see this tendency more clearly, the middle 
panel presents survival probability, p cc , as a function of 
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FIG. 3: (Color on-line) Probability to survive as cultural cen- 
ter, p cc , (middle row) and average size, (size), (bottom row) 
as a function of the distance to the nearest cultural center, 
fmin- Various lines correspond to various times t — 500 and 
t = 10 000. The average size is presented as the ratio of the 
whole system size. The top row present probability to sur- 
vive as a cultural center as a function of the position (x, y) at 
t = 500 (left map) and t = 10 000 (right map). The arrows 
in maps indicate positions of cultural centers in initial config- 
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uration. Simulation parameters: p = 0.01, 
Prepeat = 2 -6 . The system size is 100 x 100. 
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the distance to the closest cultural center, r m i n . Simi- 
larly, the bottom panel of Fig. [3] shows the average size 
of the culture as a function of r min . We can see that 
both the survival probability and average size are a non- 
monotonous function of the distance to the closest cul- 
tural center. The insertion points located very close to 
the existing cultural centers lead to maximal chances of 
surviving (p cc ) and largest size of cultures ((size)). Both 
p cc and (size) drop quickly with distance to the existing 
cultural center r m ; n and then show a slow recovery with 
the distance r m - ln . Namely, new centers either explore 
a strategy of acquiring the existing network by taking 
over a previous center, or have to explore the weaknesses 
of boundary regions to build its own network of influ- 
ence. The difficulty in building, rather than taking over 
an empire, is also reflected to the smaller size of new 
cultures emerging in distant regions, compared to new 
cultures build on deposing an existing ruler-ships. For 
larger times, both survival probability (p cc ) and average 
size ((size)) decrease for all the places (Fig. [3] middle and 



FIG. 4: (Color on-line) (A) Mean first passage time (aver- 
age survival time of a cultural center) vs distance between 
two cultural centers R in one dimensional model (filled cir- 
cles) and simplified random walk model (open circles), with 
p = 0.01. One can see good agreement for p < p rop0 at. The 
deviation grows for p > p rC pcat, but the qualitative behav- 
ior is captured in the random walk model. (B) Average size 
of cultural area in one dimension. "N" denote the numerical 
simulation and "A" denote the analytical estimate. As p r e P eat 
becomes smaller the disagreement between the simulation and 
the analysis becomes larger. 



bottom panels) since competition with existing cultural 
centers makes for a finite extinction rate. 



4- Analysis of one dimensional model 

The reported features of the present culture spreading 
model can be understood qualitatively by considering a 
simplified one dimensional model of the random walk of 
the boundary between two cultural centers. 

Suppose that there is a cultural center Co at the site 
and another cultural center C\ at the site R. We consider 
the motion of the left most site r that belongs to Co- 
Moreover, we assume that no additional cultural center 
appears, i.e. p nC w = 0. In this limit, the age a r of the 
fashion at site r which belongs to Co, and the age a r+ \ 
of the fashion at site r' = r + 1, which belongs to Ci, can 
be approximated as 



R 



1-p 



l-p 



(1) 



respectively. In Eq. ([I]), r and r' are independent dis- 
crete stochastic variables both having a probability dis- 
tribution 



P(t) — p re peat(l — Prepeat) 1 "- 



(2) 



For the sake of simplicity, the stochasticity of the fash- 
ion propagation in the preferred direction and the time 
correlation of the age are ignored, see Eq. ([I]). 

We can calculate the rate that the position r of the 
left most site that belongs to Co decreases (increases) by 
one, which happens if a r > a r +i (a r < a r +i) when the 
site r (r + 1) listen to site r + 1 (r) with probability p. 
The explicit form of the rates are given in the Appendix. 
Recall that the presented derivation is valid when p <C 
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Prcpcat, where the approximation of the age in Eq. |T]) is 
justified. 

In the left panel of Fig. |4j we compare the mean 
first passage time (MFPT) of the boundary starting at 
r = (R — l)/2 to reach one of the cultural centers, 
T(p,p repeat , i?), which gives the typical survival time of 
established cultural centers separated by distance R. The 
simplified random walk model of the boundary between 
two cultural centers agree reasonably well with the full 
one dimensional simulation results. We also checked that 
the level of agreement improves when p <C p ropoa t (data 
not shown). We can see that the survival time is longer 
for larger p rep eat and increase exponentially for large R. 

It is more intuitive to interpret these results by making 
continuous approximation for space and time and deriv- 
ing the Fokker-Planck equation for the probability P(x, t) 
that the boundary is at the position x — r — (R — l)/2 
(therefore the cultural centers are located at x ps ±i?/2) 
at time t. The resulting equations are (see the Ap- 
pendix) , 

where 
J(x,t) 

and 



dP(x,t) 
di 



dx 



J(x,t), 



(3) 



-p dX ^-P(x,t)-^ x [D(x)P(x,t)], (4) 



a probability ps p ncw R/N, where p n0 w/N is the insertion 
probability per site. Therefore, the insertion and the 
coarsening balance leads to 



T(p,Prepeat,-R c )/2 = N/(p ncw R c ). 



(8) 



The factor 1/2 on left hand side of Eq. (JsJ) comes from 
the fact that each centers are competing with two other 
centers on both sides. The comparison of the average cul- 
tural center size from the one dimensional simulation and 
the estimate based on Eq. ^ with the mean first pas- 
sage time T evaluated under continuum approximation, 
see Eq. (All, are shown in the right panel of Fig. [4] The 
agreement is satisfactory for large p re peat but becomes 
worse for smaller p rep eat ■ One of the reason for disagree- 
ment is that ^(p,Prepeat) becomes considerably large for 
small Propoat- Consequently, the potential becomes flat- 
ter and newly inserted cultural centers have higher prob- 
ability to be overwritten before the boundary reaches the 
central point (x = 0), which enhances coarsening hence 
increasing the average size of cultural area. 

In the two dimensional case, the size becomes propor- 
tional to i?^, but other parameter dependence are ex- 
pected to be qualitatively the same. The me an fi rst pas- 
sage time T is proportional to 1 /p, see Eq. (All, which 



is the time scale of the dynamics, while rapidly growing 
function of both p repe at and R. Therefore, the average 
size R^, estimated with ^ are expected to decrease with 
Pi Propoat, and Pnew, what is consistent with Fig. [2] 



U(x) 
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(5) 



5. Replaying history of Europe 
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repeat 



= with har- 
repcat )) for |z| < 



The potential U(x) has a minimum at x 
monic behavior (U(x) ps x 2 /(2£(p,p repcat )) tor \x 
£(p,Pre P ea,t)) while U (x) grows linearly with \x\ for large 
values of argument 3> £{p,p repeat))- On the other 
hand, the dependence of the diffusion coefficient D(x) 
on position is rather weak and it can be considered as 
D fa p/2. The probability p also defines the mobility, 
therefore the timescale of the random walk of the bound- 
ary is proportional to 1 /p. 

Now we can estimate the typical size of the cultural 
area or the typical distance between centers R c . Suppose 
there is a cultural center at position 0, and a new cultural 
center is inserted at a distance R. If new cultural center 
will not be inserted, the time scale when one of them will 
be overwritten by the other one can be estimated by the 
mean first passage time, T(p,p repeat , R), starting from 
the stable point x — 0. During this period, however, a 
new cultural center can be inserted between and R with 
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FIG. 5: (Color on-line) The frequency histogram, presenting 
the fraction of time for which each point constituted a cultural 
center. Parameters: p n ew = 0.004 and p rep oat = 0.015, p = 
0.01. Simulation time T — 10 s Monte Carlo steps. 

Finally, we examine the dynamics of our model on the 
European map on a 200 x 200 square lattice with a land 
mass that consists of 24 156 sites. The only constraint 
that the map give is its boundary conditions where in par- 
ticular the sea is impenetrable, and thus information can- 
not travel across seas. We use parameters where roughly 
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20 cultural centers coexists (with p = 0.01, p ncw — 0.004 
and Prepeat = 0.015; as also used in Fig. 1). Assuming 
that the fastest time it takes a rumor to cross Europe 
is about 10 years, the corresponding number of updates 
on our 200x200 will be ~ 200. In this perspective the 
subsequent snapshots in Fig. 1 correspond to 50 years, 
a timescale where changes on the European map indeed 
occurred throughout the last millennium. Fig. [5] shows 
a fraction of time for which respective sites was a cul- 
tural center, illustrating that centers tend to be more 
stable or more often established at the tips of peninsulas 
or other remote regions of the continents. That is, the 
chance to be invaded by competing fashions and cultures 
diminishes in these remote edge regions. 

Additionally, we have checked robustness of the ob- 
served patterns in Fig. [5j More precisely, we have con- 
structed frequency histograms for size dependent broad- 
cast probability p ropoat (p rcpca t ~ size/system size). This 
modification is based on a picture that new fashions come 
more frequently when the size of the culture is bigger, 
which in consequence makes a bigger culture more stable 
and convincing. The presence of such a positive feedback 
weakened the contrast between peripheries and internal 
regions. Nevertheless, distant points remained harder to 
be invaded than internal points (data not shown). Wc 
also studied the effect of short-cuts that connect two 
randomly chosen remote sites, having the possibility to 
building main roads between cities in mind. Presence of 
short-cuts further reduced the contrast between remote 
and central points (data not shown). 

It should also be pointed out that the studied model 
does not account for geographical constraints like rivers, 
mountain chains, climate and population distribution, 
which are crucial for the spread of fashions and cultures 
in real life situations. It is assumed that transmission of 
fashions are purely local and in particular that fashions 
does not travel overseas. As a consequence, the European 
simulation is more an illustration of the basic principle 
of the model than a valid simulation of available infor- 
mation highways on an ancient European landscape. 

The incorporation of mentioned constraints can signif- 
icantly change the properties of the model. Rivers and 
roads constituted information paths in pre-telegraph Eu- 
rope while mountain chains provide natural communica- 
tion barriers. Contrary to geographical landscape, the 
role of varying population density is more complex and 
less apparent. On the one hand, it is natural to imagine 
that large population density leads to larger creativity to 
start a new culture, and the number of people sharing 
the same culture also affects the ability of the culture to 
convince other people. On the other hand, it is likely that 
there is a positive feedback from cultural center to the 
local population, i.e., larger population density appears 
in places which are close to existing cultural centers. It 
would be an interesting future project to incorporate such 
an effect in the present model. 



IV. SUMMARY AND DISCUSSIONS 

We have explored a simple model for emergence and 
decline of cultural strongholds, parametrized with rigid- 
ity in local social network p, probability of emergence 
p ncw and the probability at which an existing cultural 
center broadcast fashions Repeat • 

The overall assumption of the model was the postulate 
that individuals always accept the newest (locally) avail- 
able viewpoint as their own, but obviously cannot adopt 
a viewpoint that is not available in their social surround- 
ings. Social surroundings were restricted to their 4 near- 
est neighbors on a 2D square lattice, and further biased 
with their preference for listening in the direction where 
they last obtained a useful information. As probability 
to listen in other competing directions decreases (p — > 0) 
the cultural map freezes into many small regions. 

The acceptance of a fashion according to its age only 
is an important feature of our model. In some models 
of opinion dynamics, on the other hand, a set of rules, 
which on average makes an agent accept a fashion shared 
by majority, have been adopted (e.g., voter models [I] or 
majority rules models [5]). Incorporation of such rules to 
the present model should make it more difficult for a new 
culture center to emerge and grow, though quantitative 
effect depends on the exact rules. 

It is also worth mentioning the relation between the 
present model and Axelrod model [2], which is a widely 
accepted model of formation of cultural area both by so- 
cial scientists and physicists [3]. In the Axelrod model, 
each agents has a set of opinions as a vector and interact 
with the neighboring agents according to the overlap of 
opinions: It is more likely to interact when the opinions 
are close, and when they interact the agent copies one 
of the different opinion from the neighbor to its opinion 
set. In this dynamics, the conservativeness is taken into 
account as tendency to talk to the agents that has close 
opinion, and the cultural area is formed as the agents 
align their opinions with their neighbors. In a sense, a 
culture spontaneously appears via interactions between 
agents in this model. The model can show coexistence of 
multiple cultural areas but it turned out that coexistence 
is unstable against spontaneous flipping of the opinions 
[19] , and several modifications of the model has been and 
is studied to realize stable coexistence of cultures [3] . 

On the contrary, in our model cultural area is defined 
as the area that shares the same information source. The 
random appearance of the new cultural center, which can 
be viewed as a spontaneous change of the opinion set in 
Axelrod model, is actually the important feature to keep 
the multiple cultural centers against one culture taking 
over the whole system. The key feature of our model to 
make this possible is the importance of newer informa- 
tion, which give some chance for newcomer to win against 
existing cultural centers. It can be interesting to add a 
similar feature to the Axelrod model, i.e., give some rate 
to renew opinions and value newer information more to 
verify whether multiple cultures can coexist in that case. 
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Appendix A: ID random walk model of the 
boundary between two cultural centers 

Here we derive the simplified one dimensional model of 
the random walk of the boundary between two cultural 
centers. As summarized in subsection |III 4[ two cultural 
centers Co and C\ are located at the site and at the site 
R, respectively. We analyze the motion of the left most 
site r that belongs to Co (the site r + 1 belongs to C ). 

The rate the position r decreases (increases) by one, 
W r -). r -i (W r ->. r +i), is given by the probability that a r > 
a r+ \ (a r > a. r+ i) when the site r (r + 1) listens to the 
unpreferred direction r + 1 (r) at a given time step. From 
and ^ we get (note that < r < R) 



H-2r-l 



,(l-Prcpoat) 



+ 1 



w, 



j — yr—1 



and 



P 



+ 1 



1 - 



1 



2— Prcpcat 

_ R-2r-l 

(1-Prcpcat) 

2— Propeat 



(1-Prepcat) 1- 
2— Prepeat 
R-2r-l , 



forr < 



forr > ^=2 



for r < 



(l-Propoat) 



for r > 



R-p 



2 Prcpcat 

Using above transition rates one can write the master 



equation for the probability density P Tit that the bound- 
ary is at site r at time t: 

Pr,t+1 — Pr.t = W r -l^yrP r -l t t ~ tv r ->r+l Pr,t 

Assuming that the time step and the lattice spacing 
are small, we obtain the following Fokker-Planck equa- 
tion [20] 



dP(r,t) 
dt 



d_ 

dr 

2dr 



[(W(r -4- r - 1) - W{r r + l))P(r, t)} 



2 [(W(r -> r - 1) + W(r -> r + l))P(r, t)} 



Substituting r with x = r — (R — l)/2 one gets Eqs. ^ 
to 0. 

The mean first passage time T starting from x = to 
reach x = ±R/2 is given by the closed formula [2D] 



T(p, Propeat, R) 

where 



rf.s 



R/2 



o«W / I d y 



o D(y) 



pJ7'(j/) 



dy. 



Approximating Z?(x) with p/2, see Eq. (JsJ) , one gets 
T(p,p Iepeat ,R) « -- T d S [e 2t/ ^ f'e-Wdy 

PJ-R/2 I JO 



from which Eq. ([8]) can be evaluated using Wolfram 
Mathematica. 
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